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Abstract 

In  the  light  of  the  current  and  ongoing  needs  of  DoD  and 
other  agencies  for  detection  of  subsurface  structures  and 
monitoring  of  activities  within  those  structures  there  is  a 
crucial  need  for  tools  that  allow  one  to  evaluate  the 
myriad  of  proposed  VLF/ELF  sensor  systems.  The  overall 
objective  of  this  effort  is  to  develop  analytical  tools  to 
model  and  assess  proposed  ELF/VLF  sensor  systems  for 
detection  of  underground  structures.  The  approach  taken 
is  to  a)  model  the  primary  electric  and  magnetic  fields 
generated  by  local  (e.g.  current  loops  or  electric  dipoles), 
remote  (e.g.  HAARP/HIPAS),  and  natural  sources  (e.g. 
external  noise);  b)  model  the  perturbed  fields  due  to  the 
specified  underground  structures;  c)  combine  the  primary 
and  perturbed  fields  along  with  the  external  noise  for  an 
effective  field  at  the  sensor;  d)  model  the  sensor  (the  GEM- 
2  from  Geophex,  Ltd.)  given  the  manufacturer's 
specifications  and  then  using  these  results  e)  evaluate  the 
response  of  the  sensor  system  to  the  specified 
underground  structure. 


Introduction  and  Summary 
Overview 

This  is  the  Final  Report  for  the  SBIR  effort  entitled  "ELF/VLF  Electromagnetic 
Detection  and  Characterization  of  Deeply  Buried  Targets"  sponsored  by  the  Defense 
Advanced  Research  Projects  Agency.  The  overall  objective  of  this  effort  is  to 
develop  a  analytical  tools  to  model  and  assess  proposed  ELF/VLF  sensor  systems  for 
detection  of  underground  structures.  The  approach  taken  is  to  refine  and  extend 
existing  tools  to  model: 

•  Primary  fields  generated  by 

-  Local  Sources 

-  Remote  Sources 

-  Natural  Sources 

•  Perturbed  fields  due  to  underground  structures 

•  Calculate  the  total  field  at  the  sensor  location 

•  The  response  of  the  candidate  sensor  to  the  total  fields  -  The  sensor  system 
selected  is  the  GEM-2  manufactured  by  Geophex,  Ltd. 

This  effort  is  of  significance  in  light  of  the  current  and  ongoing  needs  of  DoD  and 
other  agencies  for  detection  of  subsurface  structures  and  monitoring  of  activities 
within  those  structures.  Given  this  need  it  is  crucial  to  possess  tools  that  allow  one 
to  evaluate  the  myriad  of  proposed  VLF/ELF  sensor  systems: 

•  Claims  are  often  supported  by  sparse  technical  analysis/ tests 

•  Intercomparison  among  competing  candidates  is  impossible 

•  Optimization  is  difficult  and  costly 

The  present  effort  is  meant  to  provide  MATLAB  analysis  tools  to  allow  performance 
predictions  by  ourselves  and  the  wider  ELF/VLF  community:  DoD  agencies,  SETA 
contractors,  instrument  designers.  The  theoretical  formulas  are  given  in  Appendix 
A.  These  formulas  have  been  applied  to  a  specific  commercial  sensor  system,  the 
Geophex  GEM-2  and  the  results  are  presented.  Also  given  is  an  outline  of  a  Phase  II 
SBIR  effort  to  further  validate  an  extend  these  tools. 


-1- 


Report  Outline 

The  rest  of  this  report  is  broken  into  six  sections  listed  below: 

•  Source  Modeling 

-  Remote  Man-Made  Sources 

-  Local  Man-Made  Sources 

-  Natural  Sources 

•  Structure  Modeling 

•  Perturbed  Field  Calculations 

•  Total  Field  Calculations 

•  Sensor  Response  Modeling  Results 

•  Suggested  Phase  II  Tasks 

•  Appendix  A:  Sommerfeld  Integral  Results  for  a  three-layer  medium 
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Source  Modeling 

The  ELF/VLF  electromagnetic  sources  examined  in  this  study  are  of  three  types 

•  Remote  man-made  sources,  e.g.  HAARP/HIPAS 

•  Local  sources,  e.g.  metal-detector  loop 

•  Natural  sources,  e.g.  lighting 

This  section  examines  these  three  sources  types  and  reviews  the  results  of  the  effort 
so  far  under  this  program. 


Remote  Man-Made  Sources 

The  High  Power  Auroral  Stimulation  Observatory  (HIP AS)  and  the  High  Frequency 
Active  Auroral  Research  Program  (HAARP)  are  capable  of  exciting  plasma 
structures  in  the  ionosphere  to  act  as  elevated  electric  or  magnetic  dipoles.  These 
dipoles  can  be  then  modulated  at  specified  ELF/VLF  frequencies  to  produce 
standing /traveling  wave  electromagnetic  fields  in  the  Earth-Ionosphere  Cavity. 


HAARP-induced  auroral 
_  ionospheric  stimulation  acts 
"as  ELFA^LF  transmit  antenna 


Standing/Traveling  Wave 
ELF/VLF  Field 
in  Earth/Ionosphere  Cavity 


Ionosphere 


Figure  1:  Earth-Ionosphere  Cavity 


As  shown  in  Figure  1,  these  fields  could  then  excite  underground  structures  at  some 

distance  (thousands  of  kilometers)  from  the  source.  These  perturbed  fields  would  be 

measured  by  local  sensors.  This  section  deals  just  with  the  fields  excited  by  these 
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remote  man-made  sources. 


Perfectly  Conducting  Earth-Ionosphere  Cavity 

The  Earth-Ionosphere  Cavity  has  been  studied  by  many  researchers  utilizing  a 
variety  of  approximate  mathematical  techniques  and  any  number  of  ionospheric 
models.  Since  this  program  is  an  effort  designed  to  produce  tools  that  will  provide 
the  user  with  validated  results,  we  felt  that  a  solid  set  of  test  cases  built  on  exact 
solutions  was  required  before  developing  more  elaborate  models.  To  this  end  we 
have  examined  the  case  of  a  perfectly  conducting  spherical  cavity  excited  by  the 
following  four  dipole  sources  (Barrick,  1998): 

•  Vertical  Electric  Dipole  (VED) 

•  Horizontal  Electric  Dipole  (HED) 

•  Vertical  Magnetic  Dipole  (VMD) 

•  Horizontal  Magnetic  Dipole  (HMD) 

The  form  of  the  solution  to  the  perfectly  conducting  spherical  cavity  for  these 
sources  is  written  in  a  straightforward  way  to  yield  E/H  solutions  of  the  form 
(assuming  the  dipole/loop  is  situated  at  a  radius  b  above  the  "North  Pole"  of  a 
spherical  polar  coordinate  system)  with  the  observation  point  at  {r,6,(l))): 

cos(m[^  -  po  ])X  «« ib)F„  {kr)P„"'  (cos  6) 

«=o 


where 

b 


«n(^) 

m 


't’o 

F,(kr) 


k 


-  radius  from  center  of  Earth  to  location  of  the  dipole/loop  source 

-  coefficient  determined  by  boundary  conditions  and  dipole  height 

-  0  or  1  depending  on  field  polarization  and  dipole  orientation 

-  0  or  7i;/2  depending  on  field  polarization  and  dipole  orientation 

-  J„{kr)  or  H^{kr)  depending  on  source  height. 


(z)  -  Schelkunoff's  spherical  Bessel  function 


(z)  -  Schelkunoff's  spherical  Hankel  function 


-  wave  number  of  ELF/VLF  source 


tJ 


P^{cos9)  -  associated  Legendre  function 


The  solutions  for  the  four  dipole  sources  have  been  obtained  and  programmed  in 
MATLAB.  An  example  of  this  is  given  below  for  the  vertical  E  field  and  vertical  H 
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field  due  to  a  VED. 


Figure  2:  Fields  for  vertical  electric  dipole  source  in  a  perfectly 
conducting  Earth-Ionosphere  cavity  (Zq  «  377  Q) 

As  a  check  we  compared  these  solutions  with  an  approximate  solution  due  to  Wait 
based  on  the  Watson-Transform  (see  below).  The  comparison  is  shown  in  Figure  3. 
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The  agreement  at  50  MHz  is  good  except  in  the  near  zone  (distance  <  200  km). 
Agreement  is  not  expected  at  5  Hz  because  kr  <  1  where  Wait  s  approximation  fails. 


Figure  3:  Comparison  with  asymptotic  method  of  Wait  (VED) 
Heavy  curves  are  from  Wait,  light  curves  are  the  exact  solution 


Watson-Transform  Solutions 
The  exact  form  solution  given  above 

cos(m[0  -  ])X 

/i=0 

can  be  transformed  via  analytic  continuation  (Watson  Transform)  to  another  exact 
solution  (Wait,  1996,  Chapter  VI)  with  (roughly)  the  same  form,  but  a  different  set  of 
order  indices 

cos(m[0  -  r.  (6)4  («:r)P-  (cos(®  -  «)) 

where  the  order  values  (subscripts)  v.  take  ott  complex  values  related  to  the  modes 
of  the  Earth-Ionosphere  cavity. 

The  Watson  transform  solutions  are  useful  because  at  ELF  and  VLF  frequencies  only 
one  term  of  the  sum  is  required  for  ranges  >  1000  km  and  only  about  10  terms  are 
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required  as  close  in  as  50  km.  Typically  these  solutions  have  been  studied  and 
simplified  for  higher  frequency  regions  (VLF  and  above)  using  asymptotic 
approximations  to  the  Bessel  and  Legendre  functions  of  complex  order.  The 
commercial  package  Nlathematicu®  developed  by  Wolfram  Research,  Inc.  allows 
direct  and  exact  computation  to  any  desired  decimal  precision  of  these  functions  for 
any  complex  order  and  complex  argument.  Using  hAnthctnciticci  we  have  compared 
the  conventional  solution  with  the  Watson  transformed  solution  for  the  perfectly 
conducting  surface  case.  The  solutions  are  identical.  Using  the  Watson  transform 
solution  we  have  also  included  a  lossy  ionospheric  impedance  boundary.  These 
comparisons  are  plotted  in  Figure  4.  Note  that  the  conventional  solution  and  the 
Watson  transformed  solution  overlie  each  other. 


Figure  4:  Comparison  of  Watson  transform  and  conventional 
solution  and  the  Watson  transform  solution  for  2  lossy 
ionosphere  cases  -  all  at  10  Hz.  The  moderate  loss  relative 
impedance  is  0.0028-i0.0027;  the  large  loss  relative  impedance  is 
0.029-i0.027. 
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Ionospheric  Modeling 

In  the  discussion  above,  the  fields  were  calculated  based  on  a  perfectly  conducting  or 
impedance  boundary.  At  the  Earth  surface  this  is  a  good  approximation  for  ELF  and 
VLF,  because  of  the  high  conductivity  of  earth  and  sea  water.  This  approximation 
must  be  examined  further  for  the  ionospheric  boundary.  This  is  because  a)  the 
dielectric  properties  of  the  ionosphere  gradually  change  from  approximately  free 
space  values  below  the  ionosphere  to  something  approaching  a  perfect  conductor 
and  then  back  to  free  space  values  beyond  the  ionosphere  and  b)  the  free  space 
wavelengths  at  ELF  and  VLF  are  of  the  order  of  hundreds  to  thousands  of 
kilometers  and  so  the  ionosphere  might  appear  to  be  a  relatively  thin,  penetratable 
lossy  dielectric  layer.  This  section  examines  the  validity  of  the  impedance  boundary 
layer  approximation.  Also  examined  is  the  proper  value  of  the  effective  impedance 
and  the  effective  height  of  the  layer^ . 

For  the  cases  we  are  considering  the  dielectric  properties  of  the  ionosphere  are 
determined  by  the  electron  density  and  the  electron  collision  frequency.  The  various 
layers  of  the  ionosphere  (D,  E,  F,...)  change  as  a  function  of  time  of  day,  season  of 
year,  sun  spot  cycle  time  and  latitude /longitude. 


Figure  5:  Electron  density  distribution  and  effective  electron  collision 
frequency  for  the  ionospheric  model  used  in  this  study. 

A  full  description  would  include  all  of  these  layers  and  dependencies  and  such 


’  It  should  be  noted  that  the  Earth's  magnetic  field  has  a  significemt  effect  on  the  effective  dielectric 
properties.  This  issue  will  not  be  examined  during  this  phase  of  the  effort.  It  should  be  addressed  in 
future  efforts,  especially  in  examining  HAARP/HIPAS  measured  data. 
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models  are  available  (Budden,  1985;  Davies,  1969).  For  the  level  of  accuracy  required 
for  this  study  we  will  model  the  electron  density  as  a  single  Chapman  E-layer 
(Budden,  1985,  pp  7-14).  This  electron  density  distribution  is  given  in  Figure  5  along 
with  the  effective  electron  collision  frequency  (Budden,  1985  pl2). 


T 

R 


The  ionosphere  can  be  modeled  as  a  series  of  thin,  concentric  shells  surrounding  the 

Earth,  as  seen  in  Figure  6.  This 
7  model  can  be  formulated  in  a 

number  of  ways;  we  follow  Wait's 
formulation  (Wait,  1996,  pp.318- 
321).  The  result  of  this  analysis  is 
that  the  impedance  at  each  layer  can 
be  viewed  as  looking  "out"  of  a 
non-uniform  transmission  line. 
Thus  the  impedance  at  each  layer  is 
determined  by  the  radius  of  the 
layer,  the  thickness  of  the  layer,  the 
dielectric  properties  of  the  layer  and 
the  effective  impedance  above  the 
layer. 

In  our  case  we  will  use  layers  of 
uniform  thickness  {=h)  and  let  be 
the  radius  of  the  lowest  layer.  So  the 
impedance  at  each  layer  can  be 
written  as 

zi"’ =?((>, 

Where  the  superscript  [v]  indicates 
that  this  is  the  impedance  for  the 

vth  mode  of  the  conventional  solution  or  the  Watson  transformed  solution. 

Clearly  one  may  start  at  some  large  value  of  the  radius  beyond  the  ionosphere 
where  =  Zq  and  then  step  inwards  layer-by-layer  until  one  is  in  a  region  below 
the  ionosphere  where  the  dielectric  properties  are  again  those  of  free  space.  If  the 

calculation  results  in  a  clearly  defined  impedance  layer,  then  this  can  be  used  as  a 
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lonosphere  modeled  as  a  series  of 
layers  with  an  effective  impedance, 
Zn,  at  each  layer  interface.  As  the 
distance  above  the  Earth  goes  to 
infinity,  the  impedance  goes  to  the 
free  space  value  Zq  (~377  Q) 

Figure  6:  Layered  ionosphere  model 


boundary  condition  for  either  the  conventional  or  Watson  transformed  calculation. 
The  results  of  one  such  Watson  transformed  calculation  for  a  VED  at  100  Hz  is 
shown  in  Figure  7. 


VED  - 100  Hz  -  Zeroth  mode:  =  13.723  -  i0.219 


Figure  7:  Result  of  inward  iterated  layered  impedance  calculation 


The  calculation  is  performed  as  a  sequence  of  steps 

1.  Calculate  Vq  (Watson  transform  complex  mode  number)  for  a  perfectly 
conducting  ionosphere 

2.  Use  this  value  of  Vo  (=12.941)  to  calculate  the  effective  impedance  using 
inward  iteration  to  a  point  well  below  the  ionosphere  (in  our  case  65  km) 

3.  Use  this  updated  effective  impedance  and  an  estimated  effective  ionosphere 
height  to  recalculate  Vp  (=13.723-10.219) 

4.  Recalculate  the  the  effective  impedance  using  inward  iteration  to  65  km 

5.  Estimate  the  height  of  the  effective  impedance  layer 
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This  last  step  is  performed  by  using  the  the  equation  z|''' 
assuming  a  layer  whose  bottom(=  fcj)  is  65  km  and  whose  top(=  ^2  +^)  is  a  test 
altitude  with  an  effective  impedance  of  Z{''^ .  This  layer  is  assumed  to  have  a  free 
space  dielectric  constant  (i.e.  £2  =  £0)-  ^  values  of  heights  above  65  km  and  see 

where  this  relation  breaks  down,  i.e.  where  the  value  of  Z^''^  is  not  consistent  with 
the  free  space  dielectric  assumption.  Just  below  that  point  is  the  estimated  effective 
impedance  boundary  layer.  For  our  case  this  effective  impedance  boundary  layer  is 
at  78  km  with  an  effective  impedance  of  2.123+i7.363  Q.  Using  this  result  we  can 
now  calculate  the  field  for  this  lossy  layered  ionosphere  model  and  compare  with  a 
perfectly  conducting  ionosphere.  This  comparison  is  shown  in  Figure  8. 


VED-100  Hz 


Figure  8:  Comparison  of  the  lossy  layered  ionosphere 
with  the  perfectly  conducting  ionosphere 
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Local  Sources 

Subsurface  structures  can  be  excited  by  electromagnetic  radiation  from  local,  above 
ground  sources  such  as  current  loops  or  electric  dipoles.  At  ELF/VLF  these  sources 
can  usually  be  represented  by  one  or  more  dipole  sources:  electric  dipoles  or 


excited  by  these  local  man-made  sources  above  a  2-layer  earth  as  depicted  in  Figure  9, 
meaning  the  fields  generated  in  the  absence  of  an  underground  structure. 

Electric  or  Magnetic 
Dipole  Located  at 


Figure  9:  Two  layer  earth  model 

The  discussion  of  the  perturbed  field  due  to  the  induced  currents  on  the  subsurface 
structure  will  be  discussed  in  a  later  section. 

For  horizontal  or  vertical  (electric  or  magnetic)  dipoles  the  vertical  components  (i.e. 
along  "z-axis")  of  the  fields  Ej.  and  above  the  ground  (in  region  1)  may  be 


written  as  (see  Chew,  1995,  pp.71-79),  assuming  the  source  dipole/loop  is  at  the 
origin  of  a  cylindrical  system:  (p,0,z) 


where 


F}^  -  Vertical  component  of  the  field:  or  //,,,  the  [1]  superscript 

indicates  that  the  source  is  in  region  1 

-  Function  (may  be  constant  or  zero)  of  the  azimuthal  angle,  (j) 

K  ■  ~  K  '  (region  1)  wave  number 

p  -  4^^ 

ffiK^^iz)  -  Function  of  kp  and  k^,{= 
m  -  0  for  VED/VMD  and  1  for  HED/HMD 

ttpjPp  -  Each  may  be  ±1  depending  on  dipole  orientation/ field  component 

J,  -  Height  of  dipole  above  ground  (assumes  source  at  z  =  0) 

-  Gpneralized  reflection  coefficient  at  1-2  ground  interface 


Note:  The  integrals  described  in  this  section  and  the  Perturbed  Field  Calculations  section  are  all  given 
in  detail  with  needed  definitions  in  Appendix  A  of  this  Report. 


A  similar  relation  holds  for  the  vertical  components  of  the  fields  below  the  ground 
in  region  2  due  to  the  above  ground  source. 


0 

The  definitions  are  similar  to  above  except  that  J2  is  the  height  of  the  dipole  above 
the  2-3  interface  and  the  introduction  of  the  coefficient  Aj  related  to  the  generalized 
transmission  coefficient  through  the  1-2  interface. 


The  transverse  field  can  be  obtained  from  the  vertical  fields  via  the  following  ansatz 
(Chew,  1995,  pp.  75-76).  Observe  that  we  can  write  the  electric  and  magnetic  fields 
within  a  given  layer  as  integral  transforms: 

E(r)  =  J  dkpE{r,kp)  and  H(r)  =  j 

0  0 

where 

^r,kp)  =  E,(r,kp)  +  z^i:,kp)  and  H(r,/:^)  =  H,(r,*,)  +  z4(r,A:^) 
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Chew  shows  that 


E,' 

1 

V  —  -icofiz  X  V, 
dz 

H,. 

(0^jX£-kl 

i(0£z  X  V,  V,  — 

oz  . 

A. 

where 


V, 


The  expressions  for  E,  and  H,  may  then  be  integrated  to  obtain  the  transverse 


components  of  the  fields. 


As  an  example  of  this  type  of  calculation  consider  a  1  meter  radius  current  loop 


Figure  10:  Contour  plot  for  constant  values  of  in  dB-amp/m.  The 
source  frequency  is  3  kHz;  the  source  is  a  1  meter  radius  current  loop  (1 
amp).  The  loop  is  1  meter  above  a  two  layer  earth  with  properties  given 


in  the  text. 

(VMD)  above  the  ground  at  a  height  of  1  meter(=>  d,  =  1  meter);  assume  a  frequency 
of  3  kHz  with  ground  properties  of: 
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—  =  10  ,  <T,  =0.01mhos/m  and  -^  =  20  ,  (72  =0.02 mhos/ m 
and  the  thickness  of  the  layer  is  100  meters(=>  dj  =  101  meters). The  magnetic  field 
above  ground  can  be  written  as 


-ila^ 


U 

where  a{=l  meter)  is  the  radius  of  the  loop,  /  is  the  current  in  the  loop  (for 
simplicity  we  will  assume  1  amp),  and 

"  (K  +'‘2,tK  +  *3,) +(*., 

with 


V^o 


oj 


Carrying  out  the  integration,  we  obtain  the  contour  plot  shown  in  Figure  10  for  the 
//,j  primary  field  above  ground. 


Frequency  Dependence 

Optimizing  any  sensor  system  requires  a  knowledge  of  the  frequency  dependence  of 
the  system  and  hence  the  various  components  of  that  system.  Below  we  plot  the 
frequency  dependence  of  the  source  field  only;  in  subsequent  sections  appropriate 
plots  of  frequency  dependence  of  the  other  components  will  be  presented. 

The  plot  below  depicts  the  frequency  dependence  of  fields  from  a  VMD  source  (i.e., 
a  small  loop)  1  meter  above  a  specified  two-layer  earth  as  observed  at  two  points:  1)  a 
point  1.66  meters  displaced  horizontally  from  the  source  above  the  earth-air 
interface  and  2)  a  point  3  meters  below  the  earth-air  interface  at  a  horizontal  distance 
of  2  meters.  The  1.66  meters  is  chosen  to  reflect  the  distance  from  the  source  to 
receiver  distance  of  the  sensor  system  selected  for  analysis  (GEM-2  from  Geophex, 
see  below).  The  quantities  plotted  are  the  in-phase  (real)  and  quadrature  (imaginary) 
values  of  the  z-axis  component  of  the  B-field  in  picoTesla.  The  VMD  source  is  a  3 
amp-m^  magnetic  dipole  (consistent  with  the  GEM-2).  Again  we  use  a  two-layer 
earth  with  ground  properties  of: 

^  =  10  ,  (7,  =0.01mhos/m  and  ^  =  20  ,  cTj  =  0.02 mhos/ m 

£q  ^0 
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and  the  thickness  of  the  layer  is  100  meters. 


Figure  11:  Frequency  dependence  of  the  in-phase  (real)  and  quadrature 
(imaginary)  values  of  the  z-axis  component  of  the  B-field  for  two  points 
relative  to  a  magnetic  moment  source  of  3  amp-m^.  1)  A  point  1.66 
meters  from  the  source  above  the  earth-air  interface  and  2)  A  point  3 
meters  below  the  earth-air  interface  at  a  horizontal  distance  of  2  meters 
from  the  source 

One  surprising  result  of  these  calculations  is  that  the  below  ground  quadrature  field 
is  slightly  larger  than  the  the  above  ground  quadrature  field,  even  though  the 
distance  to  the  source  is  larger  for  the  below  ground  point  than  the  above  ground 
point  -  4.47  meters  vs.  1.66  meters.  We  do  not  have  an  explanation  for  this  at  this 
point,  although  it  might  help  in  sensor  design/ optimization. 
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Natural  Sources  -  External  Noise 

A  crucial  component  in  establishing  the  feasibility  of  a  sensor  system  are  the  signals 
always  present  in  the  natural  environment.  Lightning  strokes  around  the  earth  (100 
per  second,  on  average)  fill  the  earth-ionosphere  cavity  with  "noise"  that  has 
approximately  a  1/f  field-strength  dependence  on  frequency,  down  to  ~0.1  Hz 
(Fraser-Smith  et  al.,  1991;  Lanzerotti  et  al,  1990).  The  precise  relationship  based  on  a 
fit  to  the  experimental  data  and  referenced  to  10  Hz  is: 

where  is  the  rms  magnetic  field  average  external  noise  for  a  1  Hz  bandwidth 

in  the  atmosphere  in  picoTesla  /  -VHz  and  /  is  the  frequency  in  Hz.  This  relationship 
between  the  rms  magnetic  field  average  and  frequency  is  plotted  in  Figure  12,  below. 


Figure  12:  Plot  of  the  amplitude  of  the  rms  B-field  average 
external  noise  in  picoTesla/  VHz 
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Underground  Structure  Modeling 

The  underground  structures  will  in  most  cases  be  much  smaller  than  the 
wavelength  (even  the  shortened  effective  wavelength  due  to  the  ground  dielectric 
properties).  This  fact  allows  us  to  model  the  underground  structures  as  magnetic 
and  electric  dipoles.  In  this  section  we  will  define  and  discuss  these  dipoles  and  two 
approaches  to  calculating  these  dipoles:  NEC-3  and  an  ellipsoid  dipole  model.  For 
the  ellipsoid  dipole  model  we  give  the  values  of  the  dipole  moments  when  excited 
by  an  external  field  (electric  or  magnetic). 

Dipole  Definitions  and  the  Long  Wavelength  Approximation 

Electric  Dipole  Moment 

The  electric  dipole  moment,  p,  of  a  localized  charge  distribution,  p(r),  is  defined  as 
P  =  j(Ro  +  r)p(r)^^V 

V 

where  Rq  is  the  radius  vector  from  the  coordinate  system  origin  to  the  center  of  the 
charge  distribution  and  V  is  the  region  containing  the  charge  distribution.  This 
definition  holds  for  any  frequency. 

Assume:  l)the  structure  scale  is  much  shorter  than  the  wavelength  of  the  exciting 
radiation  (Long  Wavelength  Approximation),  2)all  dimensions  stay  comparable, 
3)material  properties  stay  finite,  and  4)the  structure  does  not  have  its  own  sources, 
then  we  can  write 
psPo^ 

So  the  dipole  moment  will  linearly  proportional  to  the  volume  of  the  subsurface 
structure;  this  will  be  seen  in  more  detail  below. 


If  the  field  is  not  static,  then  we  can  write  (using  current  conservation  and  some 
vector  calculus)  with  J  being  the  current  density  and  co  the  radian  frequency. 

p  =  ljj(r)rfV 


As  a  check  we  consider  a  localized  current  source  along  the  z-axis  from  to  ^ 
P=~Jj(*’)d  r=z—^lQ5{x)5{y)d{^^-\^dxdydz  j Q{u)=\]mi  pulse  function= 


\l  u<l 
|0  u>l 


0) 


(d(M)= Dirac -delta  function 
This  is  consistent  with  the  usual  notions  of  an  electric  dipole  moment. 
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Magnetic  Dipole  Moment 

The  magnetic  dipole  moment  is  defined  as 

“  =  ^j(^o  +  *‘)xJ(r)dV 

^  V 

and,  again,  using  the  Long  Wavelength  Approximation  and  assuming  no 
underground  sources  we  have  the  approximate  relation 

Again  the  magnetic  dipole  moment  is  linearly  proportional  to  the  volume  of  the 
subsurface  structure. 

As  a  check  we  consider  a  current  loop  in  the  x-y  plane  at  the  origin  or 

J(r)  =  ^I5{r  -  a)S(cosd) 
so 

m  =  —  |r  x^/d(r-  a)5(cos0)d^r  =  —  j  zrS(r  -  a)r^drS(cosd)sindd$d^ 

=  Im^z 

Which  is  consistent  with  the  usual  notions  for  the  magnetic  dipole  moment  of  a 
current  loop. 

NEC-3  Models 

The  most  general  approach  to  modeling  a  conducting  underground  structure  is  to 
use  a  method-of-moments  wire  grid  model.  This  approach  although  powerful  does 
require  some  care  at  low  frequencies  (NEC  was  developed  with  the  resonance  region 
in  mind).  There  are  some  technical  difficulties  for  structures  whose  feature  size  is 
<10'^  A,.  Use  of  NEC  will  be  restricted  to  structures  that  are  not  too  small,  and  where 
its  capabilities  are  required  to  adequately  describe  its  electromagnetic  response;  in 
other  cases  we  will  use  dipole  ellipsoid  models  (see  following  section). 

Ellipsoid  Dipole  Models 

Ellipsoid  Geometry 

In  this  effort  we  wish  to  model  a  variety  of  underground  structures:  rectangular 
boxes,  cylinders,  spheres.  If  the  structure  does  not  exceed  the  Long  Wavelength 
Approximation,  it  may  be  adequately  represented  as  an  oriented  ellipsoid.  An 
ellipsoid  is  defined  by  the  values  of  its  three  semiprincipal  axes:  a,,  Uj,  a^.  Below  in 
Figure  13  is  a  sketch  of  this  solid. 
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Figure  13:  Ellipsoid  Dipole  model 

47r 

The  volume  of  this  solid  is  V  =  — a^a^a^.  This  model  will  be  used  for  both  the 
electric  and  magnetic  dipole  structures. 


Induced  Dipole  Moment  Values  in  the  Static  Limit 


In  the  static  limit,  the  induced  electric  dipole  moment  on  an  ellipsoid  due  to  an 
external  electric  field  Eq  =  xEq^  +  can  be  written  as  (Stratton,  pp207-213) 


p  = 


4^1,0203 


(^slnicr  ^o)^( 


Ox 


■+y 


(^struct  ^o)^Oy 


■  +  z- 


i^srnicl  ^o)^Oz 


2e. 


2e, 


where 

^strua  '  Effective  dielectric  constant  of  the  structure 

£2  -  Dielectric  constant  of  medium  surrounding  the  structure  (see  Fig.  9) 

A^,A  ,A^  -  Functions  that  depend  on  all  three  semiprincipal  axes  -  are  ,  in 
general,  elliptic  integrals;  for  a  sphere  A^  =  A,,  =  A^= 

Clearly  p  «=  expected  from  our  Long  Wavelength  Approximation  analysis 

above. 

For  a  sphere  this  reduces  to 
P  =  ~  ^0)  E- 

Z,C2  I  ^struct  ^  ^ 

We  can  write  for  an  external  magnetic  field  Hq  =  ^ox  +  V^ox  +  Stratton, 


p258) 


m  =  - 


4^11^2  ^3 


3/to 


{^^s,ruc,-f^o)Ho, 


■+y- 


{.^^slriict  f^o)^0y 


-+x- 


{^^s,n.cr-^^o)H^ 


Oz 


Zfl2  Zfl2 


where 
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t^sinicr  ~  Effective  permeability  of  the  structure 

-  Permeability  of  medium  surrounding  the  structure  (for  our  study 
this  is  iXq) 


For  a  sphere 


m=4;za, 


3^1 

>0 


■A^o 


2)^2 


H. 


The  above  relations  for  magnetic  dipole  moments  are  only  important  when  the 
structure  has  a  permeability  different  from  this  is  not  true  in  many  cases.  But 


this  result  is  useful  in  checking  the  static  limit  (see  below)  of  the  very  important 
result  concerning  the  induced  magnetic  dipole  moment  on  a  conducting  ellipsoid. 


Induced  Magnetic  Dipole  Moment  Values  of  a  Conducting  Sphere 

Wait  (1951, 1968)  has  examined  the  frequency  dependence  of  the  induced  magnetic 


Figure  14:  Frequency  dependence  of  the  functions  used  in 
calculating  the  induced  magnetic  dipole  for  a  1  meter  radius 
aluminum  sphere 

dipole  moment  of  a  conducting  sphere  in  a  medium.  He  obtains  the  following 
result 

m  =  -2mzo^(3M  +  W)Ho 
where 
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2^,,„,„(smh7-ycosh7)+/io(sinh7-7cosh7+r^sinh7)^ 

3M+iZN-  ^^^^^^^^5ijjlj^_ycosh7)-/io(sinh7-7cosh7+7^sinh7)  “ 

with  7  =  aQ-^-ioxjji^^^, .  In  Figure  14  3M  and  3N  are  plotted  as  a  function  of 

frequency  for  a  1  meter  radius  aluminum  sphere  (c7^,=3.82xl0’mhos/m  and 
ju^,  =  1.00000065/to).  In  Figure  15  3M  and  3N  are  plotted  as  a  function  of  frequency 
for  a  1  meter  radius  iron  sphere  (cr^,  =  1.03xl0’mhos/m  and  /t^,=90/to)- 


Figure  15:  Frequency  dependence  of  the  functions  used  in 
calculating  the  induced  magnetic  dipole  for  a  1  meter  radius 
iron  sphere 


If  we  examine  the  static  limit  (i.e.  let  7^0)  and  the  permeability  of  the  surrounding 
medium  is  (Xq  we  recover  the  value  obtained  previously; 


m=4;Ki; 


'struct 


■^0 


Ho. 


_  struct  J 

For  iron  (/i^,  =  90/to)  this  static  limit  evaluates  to 

"89/io 


m=4;Kio 


92/10 


Ho=2mzo^[1.935]Ho 


=>3M+/3N=-1.935 
consistent  with  Figure  15. 

Finally  it  is  useful  to  look  at  the  frequency  dependence  of  the  scale  variation,  i.e. 
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what  happens  for  various  size  spheres  as  a  function  of  frequency.  Figure  16  shows 
3M  and  3N  for  various  radii  iron  spheres. 


frequency  [Hz] 

Figure  16:  Frequency  dependence  of  the  functions  used  in 
calculating  the  induced  magnetic  dipole  for  various  radii  iron 
spheres 

It  is  clear  from  this  plot  that  for  frequencies  >  100  Hz  the  magnetic  dipole  moment  of 
an  iron  sphere  of  radius  >  1  meter  approaches  the  asymptotic  limits  of  3M->1  and 
3iV— >0,  i-e.  m  =  -27talYi.Q. 
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Perturbed  Field  Calculations 

When  the  subsurface  structure  has  been  excited  by  electromagnetic  radiation  from 
the  above  ground  sources  (local  or  remote)  we  can  approximate  the  excited  structure 
as  a  magnetic  and/or  electric  dipole  moment.  The  value  of  this  dipole  moment  was 
discussed  in  the  previous  section.  In  this  section  we  outline  the  form  of  the 
calculations  required  to  obtain  the  above  ground  field  due  to  this  dipole.  The  field 
due  to  this  dipole  is  called  the  perturbed  field.  Figure  17  shows  a  sketch  of  the 
situation. 


Figure  17:  Induced  dipole  due  below  ground  structure  and  external  field 


The  above  ground  (Region  1)  E  and  H  fields  due  to  an  arbitrarily  oriented  dipole  in 
the  first  layer  (Region  2  of  the  sketch)  of  a  2-layer  earth  problem  can  be  written  as 
Sommerfeld  integrals.  As  before  one  need  only  formulate  the  integrals  for  and 
Hj,  and  then  obtain  the  transverse  components  via  differentiation.  These  integrals 


may  be  written  as  (see  Chew,  pp71-79) 

(0)J dk^g,[k^,k,^,k,,)j„{k^p\e-‘’^^--^^  + 


where 
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^2z 


-  Vertical  component  of  the  field:  Ejz  ^2z>  [^]  superscript 

indicates  that  the  source  is  in  region  2 

-  Function  (may  be  constant  or  zero)  of  the  azimuthal  angle,  (p 


is  the  region  2  wave  number 


P  - 

8f(^p’K'^z)  ■  Function  of  k^,  k^^ ,  and  includes  the  generalized 
transmission  coefficient  through  the  1-2  interface 


m  -  0  for  VED/VMD  and  1  for  HED/HMD 

-  Negative  z  value  of  the  bottom  of  Region  2  layer  (see  Figure  13) 

z'  -  Location  of  dipole  below  ground 

-  Generalized  reflection  coefficient  at  2-3  ground  interface 

Again  the  dipole /loop  is  situated  along  the  z-axis  of  a  (p,(/>,z)  system. 

Figure  18  gives  results  for  the  above  ground  field  due  to  a  typical  below  ground 
source  (1  Amp-m^)  -  two-layer  earth  with  ground  properties  of: 


e,  =10eo ,  cr,  =  0.0 1  mhos /m and  ,  cTj  =0.02 mhos/ m 


Figure  18:  Plot  of  the  z-axis  B-field  for  a  below  ground  VMD 
source  (1  amp-m^)  at  an  above  ground  point.  The  below  ground 
VMD  is  3  meters  below  the  surface.  The  point  above  ground  is 
1  meter  above  the  surface  and  at  a  horizontal  distance  of  3.66 
meters  from  the  below  ground  source. 
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Total  Field  Calculation 

This  section  will  review  the  relation  between  the  source  excitation  and  the  received 
signal  (open  circuit  voltage)  at  an  ideal  N-turn  loop  sensor.  The  sources  considered 
will  be  those  given  above:  remote  man-made,  local  man-made,  and  natural  sources. 
The  description  will  be  in  terms  of  a  flow  of  signal  from  source  to  receiver 
analogous  to  the  "Radar  Range  Equation". 


The  flow  of  signal  from  source  to  receiver  is  as  shown  below: 


The  first  point  to  note  is  that  the  signal  at  the  receiver  is  the  result  of  a  coherent 
sum  of  three  fields:  the  primary  field  directly  from  the  source,  the  perturbed  field 
due  to  the  underground  structure,  and  the  external  noise  due  to  various  natural 
sources.  These  have  all  been  described  earlier;  below  we  will  review  each  of  these  in 
a  summary  of  the  calculations  described  in  previous  sections. 

•  Man-Made  Source 

The  source  is  expressed  in  terms  of  an  electric  or  magnetic  dipole  moment, 
whose  location  and  orientation  over  the  earth  is  specified  by  the  user.  This 
source  may  be  local  or  remote.  Previous  sections  describe  the  calculation  of 
these  primary  fields. 

•  External  Noise 

The  fields  due  to  external  noise  will  be  coded  into  a  look-up  table. 

•  Transform  Above  Ground  Fields  to  Subsurface  Structure 

Remote  Sources:  The  above  ground  fields  in  the  vicinity  of  the  source  are 
used  along  with  the  ground  properties  to  obtain  the  below  ground  fields.  This 
is  quite  simple  since  the  wavelengths  are  long  and  hence  the  field  variation 
from  above  ground  to  below  ground  is  dependent  mostly  on  the  ground 
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properties. 

T.ncal  Sources:  The  below  ground  fields  are  calculated  in  a  straightforward  way 
using  the  Sommerfeld  integrals  described  earlier. 

•  Transform  Incident  Field  to  a  Dipole  Moment 

Model  the  underground  structure  as  an  electric  and/ or  magnetic  dipole.  This 
can  be  done  simply  by  approximating  its  shape  as  an  ellipsoid  or  by  using  a 
method-of-moments  code  such  as  NEC. 

•  Structure  Reradiation  due  to  Induced  Dipole  Moment 

Once  the  structure  is  represented  as  effective  electric  or  magnetic  dipoles,  then 
that  result  along  with  the  dipole's  location  and  the  properties  of  the  two-layer 
earth  are  used  with  the  usual  Sommerfeld  integrals  to  obtain  the  above 
ground  field  due  to  the  induced  dipole  moments  of  the  subsurface  structure. 
This  is  the  perturbed  field. 

•  Field  Sensed  at  Receive  Loop  and  Converted  to  Voltage 

At  the  receive  sensor  (assumed  to  be  an  N-turn  loop,  but  not  restricted  to  this) 
the  total  field  (a  coherent  sum  of  the  primary  field,  the  noise  field,  and  the 
perturbed  field)  excites  a  voltage  across  the  loop  that  is  then  measured. 

As  an  example  of  these  calculations  we  examine  the  following  situation  (Figure  19) 


100  meters 


(52=  0.01  mhos/m 


1  meter  radius 


83=  20£o 

03=  0.02  mhos/m 

Figure  19 


-27- 


As  in  previous  calculations  we  use  a  two-layer  earth  with  ground  properties  of: 

e,  =10eo ,  cr,  =  0.0 1  mhos /m and  gj^^OSo ,  =0.02mhos/m 

For  all  regions  =  /io  =4;rxl0”’.  The  first  layer  of  the  earth  is  100  meters  thick. 


One  meter  above  the  earth  is  a  3  amp-m^  VMD  source.  In  the  first  layer  of  the  earth 
is  an  aluminum  sphere  at  a  depth  of  3  meters  and  with  a  radius  of  1  meter.  The 
horizontal  separation  between  the  source  and  the  sphere  is  2  meters.  We  will 
examine  the  H-field  at  a  point  1.66  meters  in  horizontal  separation  from  the  source 
and  at  the  same  height  as  the  source;  the  horizontal  distance  from  this  "receive" 
point  to  the  sphere  is  3.66  meters  (2  m  +  1.66  m).  The  values  in  this  example  are 
motivated  by  the  specific  system  we  will  be  analyzing  (see  below). 

For  illustrative  purposes,  we  will  now  use  the  results  as  represented  in  Figures  11, 
14,  and  18  to  calculate  the  primary  field  (due  to  only  the  source)  and  the  secondary 
field  (due  to  the  induced  dipole  on  the  sphere)  at  the  "receive"  point.  We  will  only 
consider  the  z-axis  component  of  the  H-field  for  this  example^ .  We  will  use  a 
frequency  of  1  kHz.  All  of  these  calculations  are  done  with  the  MATLAB  code  that  is 
to  be  available  to  the  community  in  the  course  of  this  SBIR  program. 


Primary  Field 

The  primary  magnetic  field  at  the  "receive"  point  due  to  the  source  is  (at  1  kHz, 
Figure  11) 


Bprman-  =  (-65584 +i2.2400)xl0-‘" 


Recall  H. 


T> 

_ ^primary 


primary' 


^«  =  (-52-19+^0.001783)x10 


.3  Amp 
m 


Induced  Dipole  Moment 

The  magnetic  field  at  the  below  ground  sphere  due  to  the  source  is  (at  1  kHz,  Figure 


11) 


4</  =  (4696  +  i2.874)xl0''^re5/a 


=(3.737+i0.002287)xl0 


.3  Amp 
m 


The  value  of  3M+i3N  for  a  1  meter  aluminum  sphere  at  1  kHz  is  (Figure  16) 


3M+i3N=0.9961+ j0.003853 


^  The  offset  source  does  induce  horizontal  dipole  moments  in  this  case  -  these  will  be  included  in  the 
evaluation  of  the  candidate  Geophex  sensor  (below),  but  for  simplicity  will  not  be  discussed  in  this 
calculation. 


-28- 


So  the  induced  magnetic  moment  is  (with  =  1  meter) 

==-(0.02339+i0.0001048)Amp*m' 

Secondary  Field 

The  magnetic  field  at  the  above  ground  "receive"  point  due  to  a  1  Amp-m^  below 
ground  (induced)  source  is  (at  1  kHz,  Figure  18) 

5,=(397.1  +  j0.6484)xl0"'^re5/a  =»  =(0.3160+i0.00005159)xl0'^^^ 

so  for  a  m,.„^=-(0.02339+i0.0001048)Amp»m^  source,  the  secondary  field  is 

=-(9.287  +  /0.05677)xl0’'^res/fl  =>  =-(7.390+j0.04517)xl0'^^^ 

Primary  to  Secondary  Field  Ratio 

The  ratio  of  the  secondary  to  primary  magnetic  field  at  the  above  ground  "receive" 
point  for  the  previously  described  situation  at  1  kHz  is 

Jh:±y-  =  -(141 .6 +/0.8704)  x  10"® 

^primary 

If  this  calculation  is  carried  out  for  a  number  of  frequencies  we  obtain  the  results 
plotted  in  Figure  20 


Figure  20:  Plot  of  the  secondary  to  primary  z-axis  H-field  ratio 
at  the  "receive"  point  due  to  the  induced  z-axis  VMD  for  the 
configuration  given  in  Figure  19  with  an  aluminum  sphere 
(high  conductivity). 
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A  similar  plot  can  be  made  for  an  iron  sphere,  this  is  shown  in  Figure  21. 


Figure  21:  Plot  of  the  secondary  to  primary  z-axis  H-field  ratio 
at  the  "receive"  point  due  to  the  induced  z-axis  VMD  for  the 
configuration  given  in  Figure  19  with  an  iron  sphere  (high 
permeability). 
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Sensor  Response  Modeling 

The  sensor  system  to  be  evaluated  by  the  tools  developed  during  this  effort  is  the 
GEM-2  from  Geophex,  Ltd  of  Raleigh,  North  Carolina.  We  have  contacted  Geophex 
and  they  have  supplied  us  with  information  on  their  sensor  systems,  including  the 
GEM-2.  There  are  a  number  of  key  issues  that  need  to  be  addressed  prior  to  any 
sensor  system  evaluation;  these  are  listed  below  along  with  relevant  information 
for  the  GEM-2. 


Issues  Relating  to  the  Sensor  Evaluation 


Type /nature  of  the  source 

VMD  (coil  parallel  to  earth  surface) 

Type /nature  of  the  sensor 

Coil  (coplanar  with  source  coil) 

Source  cind  sensor  orientations 

Both  are  coplanar  and  parallel  to  the  earth 
surface 

Source  and  sensor  height  above  ground 

Both  1  meter  above  the  earth  surface 

Source  and  sensor  lateral  distance  from  the 
underground  structure 

Variable  -  see  results  below 

Frequency  band  of  operation 

100  Hz  -  22  kHz 

Ground  properties  and  layers/homogeneity  of 
the  earth 

Two  layer  earth  (see  Figure  22,  below  for 
material  properties) 

Underground  structure  size,  shape,  location,  ;  A  1  meter  radius  sphere  made  of  either 

and  material  properties  aluminum  or  iron.  Depth  at:  2m,  4m,  6m 

External  noise  Taken  from  Figure  12  of  this  Report  using  a  0.2 

second  measurement  period  (5  Hz) 

Sensor  internal  noise 

Assumed  negligible  for  this  analysis 

Sensitivity  of  the  sensor  system 

100  PPM  of  field  strength 

GEM-2:  System  Evaluation 

We  give  a  brief  description  of  the  GEM-2  here  with  the  overall  system  properties;  a 
more  complete  description  can  be  found  in  various  publications  given  in  the 
references  (Won,  et  al,  1996;  Keiswetter  and  Won,  1997;  Won  and  Keiswetter,  1997). 
These  publications  detail  the  specific  "bucking  coil"  technique  used  to  cancel 
(or"buck")  the  primary  field  (direct  signal  from  the  source  coil)  from  the  receiver 
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coil;  this  nulling  is  good  to  about  1  part  in  10,000  or  100  PPM  (parts-per-million)  of 
the  field  strength.  From  an  overall  system  point  of  view  the  GEM-2  can  be  viewed  as 
a  board  a  few  inches  wide  and  about  1.8  meters  in  length  with  coplanar  coils  at 
either  end,  the  centers  of  these  coils  being  separated  by  1.66  meters.  One  coil  acts  as 

z 

i 


GEM-2  at  1  meter  above  surface. 
Aligned  along  x-axis. 
Length  =  1.66  meters 


100  meters 


Conducting  Sphere 
At  a  depth  D  below  surface. 
Sphere  radius  =  1  meter 


;=  ::0  ^2 

Regions 


Sphere  Metal  Properties 

|llAi=  1.00000065po 
Gai=  3.82  X  10^  mhos/m 


For  all  regions: 
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jXFe  =  90llo 

GPe  =  1.03  X  10^  mhos/m 


Figure  22:  Sketch  of  GEM-2  sensor  evaluation  configuration 
the  source  magnetic  dipole,  the  other  acts  as  a  receiver.  The  source  coil  has  an 

effective  magnetic  moment  of  3  amp-m^  (consistent  with  the  GEM-2).  The  board 

may  be  held  so  that  the  plane  containing  the  coils  is  parallel  or  perpendicular  to  the 

earth.  The  parallel  configuration  is  called  the  horizontal  coplanar  configuration  and 

corresponds  to  a  VMD  configuration;  the  perpendicular  configuration  is  called  the 

vertical  coplanar  configuration  and  corresponds  to  an  HMD  configuration.  In  our 
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Evaluation  here  we  will  examine  the  horizontal  coplanar  (VMD)  configuration.  The 
GEM-2  can  generate  multiple  frequencies  over  a  period  of  about  0.2  seconds;  this 
permits  a  wide  variety  of  waveforms  to  be  radiated.  Calculations  are  made  by 
moving  this  GEM-2  sensor  pair  over  a  region  above  a  conducting  sphere.  The  line 
determined  by  the  sensor  pair  is  parallel  to  the  x-axis  (see  Figure  22).  We  examine  a 
number  of  configurations:  varying  frequency,  target  depth,  and  material  properties. 
In  all  cases  we  plot  in-phase  (real)  and  quadrature  (imaginary)  values  for  the  ratio  of 
the  secondary  H,  field  to  the  primary  H,  field  at  the  receiver.  The  calculations  will 
include  the  complete  effective  induced  magnetic  dipole  on  the  sphere:  both  the 
induced  VMD,  the  induced  x-axis  HMD,  and  the  induced  y-axis  HMD.  We  will  give 
the  results  as  3-D  color  plots  showing  the  secondary  H^  field  to  the  primary  H^  field 
ratio  at  the  receiver.  The  in-phase  and  quadrature  are  given  separately.  The  plots  in 
each  Figure  are  at  a  given  frequency  (100  Hz,  1  kHz,  10  kHz);  the  variation  within 
the  Figure  are  the  depth,  D,  of  the  sphere  (2  m,  4  m,  6  m)  and  the  sphere  material 
(aluminum,  iron).  The  color  intensity  is  chosen  to  threshold  at  about  100  PPM, 
consistent  with  the  sensitivity  of  the  GEM-2.  The  effect  of  average  external  noise  is 
negligible  for  this  configuration  at  these  frequencies  of  operation  as  seen  below  in 
the  calculated  noise  values  in  PPM  relative  to  the  primary  field  (using  Figure  12,  5 
Hz  BW,  and  «  0.05  Amp/m). 

Frequency  Noise  Value  Ratio 

100  Hz  1.08  PPM 

1  kHz  0.068  PPM 

10  kHz  0.0043  PPM 

Figures  23-25  examine  the  in-phase  results  and  Figures  26-28  show  the  quadrature 
results. 

The  variation  of  the  in-phase  signal  as  the  depth  is  increased  is  dramatic:  at  2  meters 
the  signal  is  "off  the  charts",  at  4  meters  the  signal  is  quite  detectable  but  greatly 
reduced,  and  at  6  meters  the  signal  is  below  the  threshold.  The  frequency  variation 
of  the  in-phase  signal  is  slight  -  the  only  noticeable  effect  is  the  slight  increase  in  the 
iron  sphere  signal  at  1  kHz  and  10  kHz  over  the  100  Hz  signal. 

For  all  but  the  2  meter  depth  iron  sphere  the  quadrature  results  are  below  threshold. 
At  100  Hz  this  signal  is  well  above  threshold  and  it  drops  to  a  marginal  value  at  10 
kHz  for  the  quadrature  signal. 
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Figure  23:  Plots  of  in>phase  secondary  to  primary  field  ratio  for 
GEM-2  sensor  evaluation  configuration  given  in  Figure  22  for 
an  aluminum  sphere  (left)  and  an  iron  sphere  (right)  at  100  Hz. 
Depth  of  sphere,  D,  varies  from  2  to  6  meters. 
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Figure  24:  Plots  of  in-phase  secondary  to  primaiy  field  ratio  for 
GEM-2  sensor  evaluation  configuration  for  an  aluminum 
sphere  (left)  and  an  iron  sphere  (right)  at  1  kHz.  Depth  of 
sphere,  D,  varies  from  2  to  6  meters. 
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10  kHz,  In-phase 
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Figure  25:  Plots  of  in>phase  secondary  to  primary  field  ratio  for 
GEM-2  sensor  evaluation  configuration  for  an  aluminum 
sphere  (left)  and  an  iron  sphere  (right)  at  10  kHz.  Depth  of 
sphere,  D,  varies  from  2  to  6  meters. 
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100  Hz,  Quadrature 

Aluminum  iron 


Figure  26:  Plots  of  quadrature  secondary  to  primary  field  ratio 
for  GEM>2  sensor  evaluation  configuration  for  an  aluminum 
sphere  (left)  and  an  iron  sphere  (right)  at  100  Hz.  Depth  of 
sphere,  D,  varies  from  2  to  6  meters. 
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1  kHz,  Quadrature 
Almninum  iron 


Figure  27:  Plots  of  quadrature  secondaiy  to  primary  field  ratio 
for  GEM-2  sensor  evaluation  configuration  for  an  aluminum 
sphere  (left)  and  an  iron  sphere  (right)  at  1  kHz.  Depth  of 
sphere,  D,  varies  from  2  to  6  meters. 
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10  kHz,  Quadrature 

Alumimim  Iron 


Figure  28:  Plots  of  quadrature  secondary  to  primary  field  ratio 
for  GEM*2  sensor  evaluation  configuration  for  an  aluminum 
sphere  (left)  and  an  iron  sphere  (right)  at  10  kHz.  Depth  of 
sphere,  D,  varies  from  2  to  6  meters. 
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Suggested  Phase  II  Tasks 

Below  we  list  some  suggested  tasks  for  a  Phase  II  effort  that  would  be  consistent  with 
further  enhancing  and  validating  the  tools  developed  during  Phase  I. 

•  Remote  Man-Made  Source  System  Evaluation  and  Validation 

Identify  a  sensor  system  that  could  be  used  with  HAARP/HIPAS  and  evaluate 
it  using  the  MATLAB  tools;  test  it  in  coordination  with  HAARP/HIPAS 
facilities. 

•  Local  Man-Made  Source  System  Evaluation,  Optimization,  and  Validation 
Using  the  MATLAB  tools,  demonstrate  techniques  for  optimizing  a  given 
sensor  system.  Extend  the  system  evaluation  of  the  GEM-2  to  other  structures, 
e.g.  pipes.  Perform  a  system  evaluation  on  the  GEM-3  from  Geophex.  Validate 
the  results  from  these  evaluations  and  optimization  using  controlled 
experiments.  Examine  the  feasibility  and  sensitivity  of  Geophex-derived 

sensors  for  aircraft  surveillance. 
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